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The basic equations for the description of structure formation in cosmology are Einstein's 
laws for gravitationally interacting systems or, with restrictions, Newton's law of gravity. 
Taking the Universe to be filled with a self-gravitating pressure-less fluid ('dust'), one has to 
deal with a nonlinear system of differential equations. To simplify this system of differential 
equations one commonly assumes special symmetries. Then, for each set of symmetries, we 
, obtain a class of solutions. 

I For simplicity and for some philosophical reasons, most cosmological models studied today 

■ are those based on the assumption of homogeneity and isotropy. Observationally one can find 
. evidence that supports these assumptions on very large scales, the strongest being the almost 
I isotropy of the Cosmic Microwave Background radiation after assigning the whole dipole to 

■ our proper motion relative to this background. However, on small and on intermediate scales 
up to several hundreds of Mpcs, there are strong deviations from homogeneity and isotropy H . 
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Here the problem arises how to relate the observations with the homogeneous and isotropic 
models. The usual proposal for solving this problem is to assume that Friedmann-Lemaitre 
models describe the mean observables. Such mean values may be identified with spatial 
^ I averages. For Newtonian fiuid dynamics the averaging procedure has been discussed in detail 

' in • The key difference between the averaged models and the standard model can be related 

to the fact that the spatial average of a tensor field A over a domain comoving with the fiuid 
^ , (a "Lagrangian domain" ) does not commute with the time evolution of that tensor field: 

dt{A)v - {dtA)v = {Ae)v - {e)v{A)v ■ (1) 

{A)t> denotes the Euclidean spatial average of A over a domain V, and 9 the local expansion 
rate. 

By averaging Raychaudhuri's equation for a dust matter model and using the commutation 
rule (g) one obtains a differential equation for the averaged expansion rate. This equation 
can also be written in a form similar to the standard Friedmann equation, but now with the 
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domain dependent scale factor a© = ; Vxi = \'D\, featuring an additional 'backreaction' 



term Q: 



3— + 4TTG{g)T> -A = Q with Q := l{{9^)v - {0)1) + 2{lo^ - a')v ■ (2) 



a denotes the rate of shear and uj the rate of rotation of an infinitesimal fiuid element. 

As soon as inhomogeneities are present in the domain, the rate of shear, the rate of rotation 
and the expansion-rate are nonzero. Hence, Q ^ and the domain-dependent scale factor 
ax) will behave in a different way compared with the scale factor of a homogeneous-isotropic 
Friedmann cosmology. 
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1.2 Spatially compact cosmologies without boundary 

The backreaction term can also be written in the following form after the application of GauB' 
theorem: 

Q = ^f (u(V-u)-(u.V)u).dS-|(-^ / u-ds)', (3) 

Vv JdV O \Vt> JdV J 

with the surface dV bounding the domain T) and the surface-element dS; u denotes the 
peculiar velocity field. 

If we choose the compact domain T> to be the whole universe, then backreaction vanishes 
for all models with closed 3-spaces (compact without boundary). For example all toroidal 
models with local inhomogeneities can be treated globally like Friedmann models. The ho- 
mogeneous density in the Friedmann models is then simply the density averaged over the 
whole space. From the same argument we conclude that no backreaction is present in A^- 
body simulations on the periodicity scale. However, this can be accomplished only within the 
framework of Newtonian theory. In General Relativity the situation is more involved due to 
the presence of an averaged contribution from the Ricci curvature, and the line of arguments 
above is not conclusive in that case Q. 

1.3 Effects of backreaction on intermediate scales 

On intermediate scales the Newtonian 'dust' approximation is well-established in models of 
structure formation. With the inhomogeneous Friedmann equation (^) we have a tool to 
quantify the time evolution of the scale factor of spatial domains in the Universe and to relate 
them to the global background expansion. Unfortunately, the dynamical evolution of the 
backreaction term is not known. Therefore, we use the Eulerian linear approximation and the 
'Zel'dovich approximation' to estimate the effect of the backreaction term. Such a calculation 
is only well-defined, after all what has been said, if we assume that the Universe can be 
approximated by spatially flat space sections and that the inhomogeneities are subjected 
to periodic boundary conditions on some large scale. In this way the calculation of the 
backreaction effect on scales below the periodicity scale can be investigated. We want to 
emphasize that no conclusion about the global value of backreaction is possible, because it 
vanishes by assumption. This, of course, restricts the generality and shows the need to go 
to a general relativistic investigation. Russ et al. |^] have attempted this in a recent work. 
However, their assumptions to start with are too restrictive, so that a vanishing backreaction 
already follows from their basic equations (see |3| Q). 

1.3.1 Backreaction in Eulerian linear theory 

We use the linear approximation theory to calculate the time evolution of the velocity field. 
For an Einstein-de-Sitter background and in the initial stages of structure formation the 
backreaction term Q decreases 

/ + \ 2/3 

Q^ = a-iQo with a= (^-j . (4) 

The dimensionless contribution to the expansion may be estimated by Q/ {^t^Gqh)^ quanti- 
fying the impact of the backreaction Q in comparison to the background density qh- This 
dimensionless backreaction grows proportionally to the variance of the density contrast field 
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'^r.m.s. ^ ^ • For this type of linearization we have to assume that the deformation of the 
comoving volume is negligible. We can relax this assumption if we work in the Lagrangian 
picture. That is why the 'Zel'dovich approximation' is a better tool in this context. 

1.3.2 Backreaction in the 'Zel'dovich approximation' 

Zel'dovich's approximation is a subcase of solutions to a Lagrangian first-order perturbation 
approach for the (dimensionless) trajectories of the fluid elements we have: 

/^(X, t) = a(t)(X + eWVoV'lX)) . (5) 

Here, VoV'(X) denotes the initial displacement field, ^(t) is a universal function of time and 
X is the initial position of the fluid particles. The 'Zel'dovich approximation' is known to 
give a good picture for structure formation also in the nonlinear regime until shell-crossing 
takes place. 

It is useful to calculate the scaled backreaction term a|,Q in the 'Zel'dovich approxima- 
tion', which can be written in a compact form [Q]: 

a%Q = a^i^ (2(//oK-|(/o)|,„ + 

with {IIo)voj {IIIo)vo being the first, the second and the third scalar invariants of the 

tensor field ip^ij averaged over the initial domain I^(to); '|' denotes the derivation with respect 
to Lagrangian coordinates. 

Using the approximation (P) for the backreaction Q we solve the differential equation (|2|) 
for aj) numerically. The results for various initial displacement fields are shown in Fig. [l| for an 
Einstein-de-Sitter background model with a = {t/tof'^'^ . We have chosen {Io)vo = in both 
plots, corresponding to a mean over-density of zero in the domain Vq. If no backreaction is 
present, such domains should follow the expansion of the background model. However, with 
(/Io)x>o 7^ and {IIIq)t>q 7^ we have Q / in general. The accelerated expansion visible 
in the left plot and the accelerated collapse in the right plot is triggered only by the small 
deviations in the initial displacement field as described by {IIq)dq and {IIIq)x>o- 

Summarizing, we have shown that the backreaction significantly influences the dynamics of 
Lagrangian domains in the 'Zel'dovich approximation'. Backreaction can invoke a dynamics 
that ressembles cosmologies with a cosmological constant. Indeed every homogeneous and 
isotropic A-model can be approximated with a backreaction model for suitably chosen initial 
conditions. Considering Gaussian random fields we will estimate the effect for generic initial 
displacement fields Q. Backreaction may, on the other hand, provide us with a new class of 
collapse models that outperform the standard 'top-hat' model (which excludes backreaction 
by its restriction to spherical symmetry). Such models may deepen our understanding of the 
dynamics on cluster scales. Finally, backreaction may substantiate the problem of missing 
mass: for domains dominated by shear fluctuations, backreaction can act as a "dynamical 
dark matter" component. 
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Figure 1: In the left panel the time evolution of the scale factor a-p is shown for {Iq)t>q = 0, 
(/Io)p^ = 10-6, and {IIIo)vo = 10"^ (dotted), {IIIo)vo = 10"^ (short dashed), {IIIo)vo = 
10~^ (long dashed). In the right panel we used {Io)vo = 0, {IIo)vq = — IQ-^, and {IIIo)x>o = 
-10-^ (dotted), {IIIo)vo = -10"'^ (short dashed), (//Io)do = -10"^ (long dashed). The 
solid line corresponds to an Einstein-de~Sitter universe, i.e. Q = 0. 
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